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Abstract. We briefly discuss the method of population synthesis to calculate theoretical delay 
time distributions of type la supernova progenitors. We also compare the results of the different 
research groups and conclude that although one of the main differences in the results for single 
degenerate progenitors is the retention efficiency with which accreted hydrogen is added to the 
white dwarf core, this cannot explain all the differences. 
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1. Introduction 

Given the uncertainties in the theoretical derivation of what the progenitors of type la 
supernovae (SNIa) are, statistical methods can in principle be useful. Viable progenitor 
scenarios must not only potentially produce a SNIa (i.e. an explosion that looks like the 
observed SNIa), they must also occur often enough to explain the observed/inferred 
SNIa rate. In order to do so, an estimate needs to be made about the occurrence 
frequency of di fferent types of binaries. This can be done using population synthesis 
techniques (e.g. IPostnov and Yungelsonl |20Q6|) as has been done by several groups for 



the different SNIa progenitor scenarios (e.g. j Yunge" 



Han an d Podsiadlowski 200J iForster et al.l 12006 



son and Livio 1998; Yungelson 2005; 
Ruiter et al. 2009; Mennekens et all 



2010: IWang et al.ll2009l " |20M . In the following we will briefly discuss some of the impor- 



tant ingredients of population synthesis calculations, uncertainties and possible mitigat- 
ing efforts before showing a comparison of results of different groups and an attempt to 
understand the differences. 



2. Population synthesis 

The basics concept of population synthesis is simple: for a (large) sample of initial bi- 
naries (i.e. two Zero- Age Main Sequence stars orbiting each other with a certain orbital 
period) the subsequent evolution is determined. Folding those with assumed initial pa- 
rameter (masses, periods) distributions and if needed star formation histories, any binary 
population can be synthesised. 

Population synthesis calculations typically evolve so many binaries that using a stellar 
evolution code that solves the structure equations explicitly is not feasible. Even if that 
would be the case, there are many binary evolution scenarios that cannot be computed 
with such codes anyway. Therefore people use rapid fits to evolutionary calculations or 
interpolations in pre-computed grids. This introduces uncertainties, although these are 
typically small compared to other uncertainties. 

More problematic is determining the outcome of binary interactions. If one star in 
the binary fills its Roche lobe and starts to transfer matter to the companion, one first 
has to determine if this will lead to runaway mass transfer and the (likely?) start of a 
"common-envelope" phase. If not, the mass transfer is deemed "stable", but it still needs 
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Figure 1. Period - mass distribution of the observed double white dwarfs compared to the 
results of the population synthesis of Toonen et al. in prep. The simulation matches the obser- 
vations reasonably well. 



to be determined how much of the transferred mass is accreted by the companion (likely 
depending on the speed at which the mass is transferred) and how the remaining mass 
leaves the system (i.e. how much angular momentum is lost along with the mass). That 
will have a large impact on the further evolution of the system. 

In the context of SNIa, the two main progenitor scenarios are carbon-oxygen white 
dw arfs that either m erge with another carbon-oxygen white dwarfs (" double degenerate" 
e.g. lWebbinkl dl9 84^ or are accretin g from a non-degenerate companion ("single degen- 
erate e.g. IWhelan and Ibenl (|l973h ). The most important uncertainties for population 
synthesis of these populations are the outcomes of mass transfer and for the SD progen- 
itors in addition the question which white dwarf — companion star configurations lead 
to sufficient accretion to reach ignition conditions by approaching the Chandrasekhar 
mass. These uncertainties are sufficiently large that we simply assume or hope that the 
prescriptions we use about these mass transfer phases in our population synthesis codes 
are right. Instead we have to use available observational data to calibrated the models. 

2.1. Comparison with observations 

Comparison of the simulated populations with the observed population of (carbon- 
oxygen) white dwarfs with non-degenerate companions has not be done yet to our 
knowledge. One of the reasons is that the observed population is severely biased by 
observational selection effects that need careful modelling before the observations can 
be compared to the population synthesis calculations. We are planning to do that in a 
forthcoming paper. The population of close double white dwarfs is much more homoge- 
neous and we have compared the results of our population synthesi s calculations with the 
observations to constrain the uncertain mass transfer phases (e.g. iNelemans et al.l 12001 
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Figure 2. Rescaled DTDs for the DD scenario for the different groups. 

Toonen et al. in prep). In Fig. [1] we show this comparison for our newest simulation of 
the double white dwarf population (Toonen et al. in prep.). 

There is a limit to the accuracy of this comparison when considering DD progenitor 
scenarios, because non of the observed systems is likely a double carbon-oxygen white 
dwarf with sufficient mass. However, the agreement with the other double white dwarfs 
does give some confidence in the simulation results. 



3. Delay time distributions of different groups 

Last year for the Lorentz center workshop " Observational signatures of type la super- 
nova progenitors" , we asked a number of different research groups to give us their data 
on the type la supernova progenitors so that we could do a comparison of the results 
of the different groups. All groups kindly accepted the invitation and were very help- 
ful in clarifying any queries regarding their methods and normalisation. We decided to 
make the comparison of the so-called "Delay Time Distribution" (DTD), which shows 
the SNIa rate per unit mass as a function of time since a hypothetical instantaneous 
starburst. For binary population synthesis calculations this is easy to compute, the only 
issue hampering direct comparison between groups is the normalisation. This depends on 
the assumed Initial Mass Function, the percentage single stars and the initial distribution 
of mass ratios and orbital perio ds. The groups involved are 

• Lev Yungelson's work (e.g. lYungelsonll2Q10l 12005c I Yungelson and L ivio 2000|) 



• The Yunnan g roup (Han, Wang, etc), see lWang et al.l (|2010l ): lHanl (|2008l ): lHan and Podsiadlowski 
(20041): iHa^l (|l998h 

• T he StarTrack code, Belczynski/Ruiter work (e.g. Ruiter et al.ll2009l : lBelczvnski et al. 

1qq3) 



The Brussels group (see I Mennekens et al' l2Q10h 
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Figure 3. Rescaled DTDs for the SD scenario. 



• The Utrecht group, Claeys, Pols (e.g. Claeys, this proceeding s) 

• Our SeBa code (Toonen et al. (in prep), [Nelemans et al.ll200ll ; |Portegies Zwart and Verbunt 



1996) 



For the comparison we res caled the resu l ts (if needed) to the following assumptions: 
flat mass ratio distribution, a lKroupa et al. 1 (Il993h IMF, flat mass ratio distribution, flat 
distribution in logP or b and 50% binaries. The conversion factors for the different groups 
are given in table [TJ 

In figures [2] and [3] we show the (binned) DTDs for the different groups, for the DD and 
SD scenarios respectively. For the latter w e only considered systems with hydrogen rich 
donors, excludin g helium rich dono rs (e.g. IWang et al.l 12009). We also plot the inferred 
data points from lMaoz et al.l ( 2010f ). Most striking is that the double degenerate channel 
shows roughly the same shape in all the calculations, although with differences in the 
actual rates, but that the single degenerate DTDs are all over the place. 

Because it is rather hard for the eye to integrate the area in a log - linear plot we show 
in table [T]the integrated rates over a Hubble time (i.e. in total number of SNIa per M in 
a Hubble time), together with the inferred rate from M aoz et aD (|2010l ). It is clear that 
all population synthesis calculations (with this normalisation!) produce too little SNIa. 



4. Can we understand the differences? 

We investigated the most likely cause of the enormous discrepancy between the groups 
for the SD scenario: the actual conditions in which an accreting white dwarf can grow in 
mass and thus explode as a SNIa (see Bours et al. in prep for more details). The different 
groups use a different approach: some use so-called retention efficiencies (the fraction of 
accreted hydrogen that in the end is added to the core of the white dwarf as carbon and 
oxygen), some use progenitor "islands" in parameter space (of orbital period, white dwarf 
and companion mass). We first calculated the SD DTD by using the islands as calculated 
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Table 1. Normalisation factors and integrated SNIa rates in 10 _4 M© _1 for the different groups 
and the diff e rent c hannels. For comparison the integrated rate inferred from the observations by 
iMaoz et"aTI (|201Q[ ) is also given. 



Group 


DD 


SD 


factor 


Yungelson 


2.4 


0.006 


0.8 


Wang/Han 


4.4 


2.8 


0.77 


Ruiter 


5.7 


0.17 


1? 


Mennekens 


2.2 


3.7 


0.62 


Claeys 


1.9 


0.13 


1 


Toonen/Bours 


1.9 


1.1 


1 


Observed 


23 
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Figure 4. Retention efficiencies (i.e. efficiency with with accreted hydrogen is tuned into 
carbon and oxygen used by the different groups. 



bv lHachisu et al. (2008), the same approach as taken bv lMennekens et all (|2010h . We also 
im plemented retention efficiencies that are the basi s of these is l ands, ones that are used 
bv IRuiter etHI (|2QQ9h and ones that are used bv lYungelsonl (I20T0I ) (see Fig. ED) . The 
latter two are based on calculations by IPrialnik and Kovetzl ( 1995 ): [~ 



); llben and Tutukov 

(1996). The resulting DTDs of these calculations are shown in Fig. As can be seen, 
the different retention efficiencies do cause a large difference in the integrated rates (with 
no systems left with the "Yungelson" efficiencies), as is the case in the comparison in 
the previous section. However, the differences in the shape of the DTDs as found in 
the previous section are not reproduces in our modelling. This suggest that apart from 
the retention efficiencies other differences in the codes, most likely in the treatment of 
the mass transfer phases contribute significantly to the uncertainties in the theoretical 
DTDs. 
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Figure 5. DTDs simulating the different app roaches: "is l ands" for assuming the parameters 
that lead to SD SNIa to lie in the islands of Hach isu et all (12008T) . The others are for th e reten - 
tion efficienc i es sh own in Fig. [4] used bv lYungelsonl (2010, 2005); Yungelson and Livid ([2000), 
Ruit er et al.1 (2009) and one that should reproduce the islands. 

5. Conclusions 

We have shown that for the theoretical DTDs for the DD scenario there are (limited) 
observational tests that can be made to constrain the results. Interestingly, the DTDs of 
different research groups for the DD scenario agree reasonably well. For the SD scenario 
the DTDs differ wildly. This is partly due to the extreme uncertainty of the retention 
efficiency of accreted hydrogen on carbon-oxygen white dwarfs. However, a test in which 
this is varied while keeping the rest of the population synthesis code fixed, shows that 
this alone cannot explain the differences in the SD DTDs between the different research 
groups. Before the SD DTDs can be used to make any statements about the likelihood of 
this scenario, the retention efficiencies have to be determined better and the progenitor 
population has to be compared in detail to the local observed sample. 
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